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The 2D lattice gas model with competing short and long range interactions is used for calculation 
of the incoherent charge transport in the classical strongly-correlated charge segregated polaronic 
state. We show, by means of Monte-Carlo simulations, that at high temperature the transport 
is dominated by hopping of the dissociated correlated polarons, where their mobility is inversely 
proportional to the temperature, /i oc T _1 . At temperatures below the clustering transition tem- 
perature the bipolaron transport becomes dominant. The energy barrier for the bipolaron hopping 
is determined by Coulomb effects and is found to be lower than the barrier for the single-polaron 
hopping. This leads to drastically different temperature dependencies of mobilities for polarons and 
bipolarons at low temperatures. 



I. INTRODUCTION 

The physics of small and large polarons is increasingly 
one of the most intriguing branches of contemporary solid 
state physics. Polarons are often discussed in connec- 
tion to the transport and optical properties of doped 
semiconductors^. The kinetic properties of polaronic 
systems are well studied in the limit of small density of 
polarons^. However, the situation becomes more compli- 
cated in the case of heavily doped semiconductors and in- 
sulators. The correlations between polarons may lead to 
the formation of a homogeneous gas of bipolarons or the 
formation of a charge ordered^ or phase separated stated. 
Indeed charge segregation of polarons has recently been 
discussed in literature^. The kinetic properties of a di- 
lute bipolaron gas are known in the band limit 7 - and in 
the limit of hopping conductivity^. On the other hand 
the kinetic properties of charge ordered and especially 
phase separated polarons and bipolarons remain an open 
and challenging problem. 

The generally accepted mechanism of charge transport 
in doped semiconductors at low temperature is variable 
range hopping (VRH) of localized electrons in the im- 
purity band^. Inclusion of the long range Coulomb in- 
teraction leads to the opening of the Efros-Shkolvskii 
Coulomb gap and a modified temperature dependence 
of the conductivity, compared to the VRH models. In- 
creasing the complexity further, there have been in- 
creasing debates about the role of multi-electron hop- 
ping conductivity. Tielsen and Schreiber have per- 
formed extensive calculations and came to conclusion 
that the many-electron hopping processes become dom- 
inant at low temperatures^, which was also supported 
by simulations^. Later, these results were challenged 
by Tsigankov and Efrosi^ who argued that a different 
choice of hopping rates leads to a substantial suppres- 
sion of the two-particle contribution. More recently^ 
it was shown that if one assumes more general hop- 
ping rates, the many-particle contribution to the con- 
ductivity may remain important at low temperatures. 
Thus the many-particle contribution to the conductiv- 



ity of doped semiconductors and insulators, particularly 
at low temperatures remains controversial and limits our 
understanding of the transport not only in doped polar 
semiconductors, but also doped cuprates (including high- 
temperature superconductors)^, doped manganitesS and 
many other materials^ where the electron-phonon inter- 
action is intermediate or strong. 

Recently we have shown that doping a semiconductor 
with charged carriers leads to the segregation of polarons 
on short length scaled. We have shown that polaron clus- 
tering with an even number of particles is more favorable 
than with an odd number of particles, because of the par- 
ticular symmetry of the short range potential and topol- 
ogy of the lattice. This suggests that we can consider 
these clusters to be effectively composed from bipolarons 
and it is natural to expect that the low temperature con- 
ductivity will be determined by bipolaron hopping be- 
cause the formation of clusters with an odd number of 
polarons are energetically less favorable. In our model 
calculations^ the distance between polarons within the 
bipolaron is small (inter-site bipolarons) so the model 
calculations of electronic transport are less sensitive to 
the choice of the hopping probability^ than in the case 
of the standard Coulomb gas. 

Here we present extensive numerical simulations of 
polaroni 7 - and bipolaron hopping conductivity in the 
charge segregated state which appears as a result of dop- 
ing the strongly-correlated insulator. We show that at 
low concentration and low temperature bipolaron trans- 
port becomes dominant whereby the addition of a short 
range attraction caused by the lattice deformation to the 
standard Coulomb gas Hamiltonian unambiguously leads 
to the dominant role of the two-particle conductivity at 
low temperatures. 

II. THE MODEL 

To perform the simulations we chose the lattice gas 
model on a 2D square lattice with competing long-range 
Coulomb and short range anisotropic Jahn- Teller (JT)- 
like interactions, because it naturally describes the for- 
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mation of bipolarons and bipolaronic clusters. The model 
in terms of the pseudospin operators, S*, can be written 



limit of the vanishing electric field to: 
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H = £(V JT (i - j) S?S? + V C (i - j) QiQj), (1) 



where S? = 1, S? = — 1 represent a single particle in ei- 
ther of the states of the electronic JT doublet at the site 
i and S? = represents an unoccupied site. Q\ — (S?) 2 
is the on-site charge and Vc (m) = 1 Jm is the dimen- 
sionless Coulomb potential. Vjt (i — j) represents the 
anisotropic short range JT interaction^^. This Hamil- 
tonian was used to demonstrate the stability of a variety 
of different textures which might appear upon doping of 
strongly correlated 2D insulators. It was shown^^ that 
the short range attractive potential alone leads to global 
phase separation which is frustrated by the long-range 
Coulomb potentia l 18 ' 19 . The system has a glassy ground 
state with many minima in the free energy. 



A. Single polaron hopping 

To analyze the charge transport we use standard vari- 
able range hopping transport theory 2 ^ The applicabil- 
ity of this theory is justified by the glassy ground state 
demonstrated by previous simulations^^. The single 
particle hopping rate between an occupied site i with 
pseudospin a to an empty site j with the resulting pseu- 
dospin P is given by>2£ 



Aiaj/3 > 
Aiaj/3 < 
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where A\ a ^ = AH (ia — >jy8) — Exij is the dimension- 
less many body energy difference between the final and 
initial states in the presence of an external electric field, 
E, applied along the a;- axis. T is the temperature and 
7ocxp (— 2a sp nj) is the transition probability, r,j = |j — i| 
and a~p the spatial extent of the single polaron wavefunc- 
tion. For a given system configuration, {S^}, the total 
dimensionless instantaneous current is the sum over all 
possible hops: 
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where = j x —ix, with the first sum running over all oc- 
cupied states with S? ^ and the second over all empty 
sites with Si = with both possible final pseudospin 
states, j3 = ±1. Since only the hops with Ai Q j^ > con- 
tribute to the net current the formula simplifies in the 
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where N e g is the total number of particles/holes 2 ^ 
and /i sp the dimensionless single polaron mobility. 
The second sum runs only over the final states with 
AH (ia -> j/3) > 0. 



B. Bipolaron hopping 

Bipolaron hopping is a second order process^ where a 
bipolaron (ia, i'a') hops to (j/3, j'/3') via the virtual disso- 
ciated states i'a)' or (ia, j'/3'). Tsigankov and Efrosi^ 
argued that the hopping probability is suppressed when 
the distance between the particles within a hopping pair 
is large, which reduces the biparticle conductivity by or- 
ders of magnitude. This is not an issue in our case since 
the size of our bipolarons is restricted to one lattice con- 
stant. 

The matrix element for the bipolaron hop is given by: 
< ia,i'q'|T>|j/3,iV X j/3, i'a'\V\jP, j'/3' > ft ~ 

where the operator V is determined by Eq. (4.13) of 
Refill and has both diagonal and off-diagonal elements 
in phonon occupation numbers. Each matrix element in 
([5]) contains a factor exp (— a sp R), which depends on the 
distance of the hop. The dependence of A on the hop 
length is therefore exp [— a sp (nj + H'j')]- 

Another question which should be addressed is the fre- 
quency of the hopping events which is determined by the 
pre-exponential factor of the matrix elements. Here we 
are interested in the case when the single polaron hopping 
is allowed so the binding energy of the pair, A, is small, 
A << cj, where uj is the phonon frequency. We can there- 
fore assume that E\ a ^ a i — E- } ^^i a i ~ —uj and the pre- 
exponential factor of A is t 2 /uj, where t is the frequency 
for the single polaron jump. Since we are interested in 
the single phonon processes the second matrix element 
in ((5|) must be diagonal in the phonon occupation num- 
bers. The single polaron jumps occur with the frequency 
t ~ uj so up to a numerical factor of the order of 1 the pre- 
exponential factor of A is ~ u) , similar as in the single po- 
laron hop case. The bipolaron-hopping-rate temperature 
independent pre- factor is therefore ^ 3p e _2Qsp ( rij+, ' i 'j') j 
with 7g P w 70, and the bipolaron hopping rate is given 
by: 
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where (ia, i'a') denotes the nearest neighbors and 

&(u>,v a >)<iwfi>) =AH((ia, i'a') -> (j/3, - (?) 
-E(x )i +x i >y) 

is the many body energy difference between the final and 

I 



initial states. 



The total bipolaron instantaneous current for a partic- 
ular system configuration in the limit of vanishing electric 
field is similar to the single polaron case, 
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(io,i'a') <j/3j'/3'> 



The first sum runs over all occupied bipolaron and the 
second over all possible unoccupied bipolaron states with 
AH ((io-jiW) -4- (j/3,j'/3')) > 0. For the sake of com- 
parison with the single-polaron mobility Eq.® also de- 
fines the effective dimensionless bipolaron mobility ^ p 
per single particle/hole. 



III. RESULTS AND DISCUSSION 

We calculate the polaron and bipolaron mobilities at 
any given T by averaging equations {3]) and {8} over 
Markov chains of the system configurations obtained by 
means of MC simulations. Contrary to ref. [l3| a re- 
alistic MC dynamics, which includes correlated biparti- 
cle hops in the external field, is not required. Instead, 
unrestricted-length single-particle hops were used to im- 
prove convergence. The details of the used MC algorithm 
were described in ref 

For the present calculations we used periodic bound- 
ary conditions on a system size L x L = 30 x 30 after 
carefully checking for any size effects. The shape of the 
Vjt(i) was taken to be nonzero only for i = 1 and was 
therefore specified by a single parameter Vjt(1, 0) = Vjt- 
To minimize further the computation time we limit the 
maximum hopping distance in equations Q and (|8"j). Af- 
ter a careful check for convergence we set a sp = 0.5 and 
selected the maximum hopping distance, rc = 8. In Fig. 
[T] we plot calculated temperature dependence of mobili- 
ties at different dopings n = N/L 2 , where N is the total 
number of polarons in the system and n < 0.5 represents 
particle and n > 0.5 hole doping regions. 

At high T the second exponent in the sums Q and © 
shows a negligible temperature dependence so the dom- 
inant contribution to temperature dependence of both 



100 



10 



0.1 

100 



10 



0.1 
0.01 
1E-3 



. ■ ■■■^s. B--. a) 

* **** "^^^v. 

■ / ■ 


b) 


7 If • 0.02 
// —^0.05 
V —-—0.10 
0.20 


jff ——0.98 
*ff — — 0.95 

y — 0.90 

0.80 ' 


■ — •— 0.30 ' '■ 
*v 0.40 
^^■^ . 0.50 


■ — •— 0.70^ 
v 0.60 

— •— 0.50 ■ 


v f 


/J / ^S^V 1 " ' ** 


C): 


d) 



0.1 



4 0.1 



T 



Figure 1: Temperature dependence of the single polaron mo- 
bility, £i sp , (a) and (b) and the effective bipolaron mobility, 
/i bp , (c) and (d) at Vjt = — 1 at different doping levels n. 
In the insets the corresponding hopping processes are shown 
schematically. The solid and dashed black lines show 1/T and 
I/T 2 slopes, respectively. 



mobilities comes from the 1/T prefactor. When the tem- 
perature is decreased we first observe the suppression of 
mobilities due to the orbital (pseudospin) ordering£ii£ 
in the hole doping region, n > 0.5, which is most pro- 
nounced in the single polaron mobility. The suppression 
is due to the high energy cost of the hops with the pseu- 
dospin flip in the orbital ordered system which makes half 
of the final states almost inaccessible. 

With further decrease of temperature, the behavior of 
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the mobility at low particle and low hole doping becomes 
almost identical: The single polaron mobility shows a 
maximum at any doping at the temperature T^, which 
is doping dependent. In contrast, the maximum is absent 
for the bipolaron mobility at low particle/hole doping. 
The bipolaron mobility even starts to increase faster than 
1 /T at low temperatures exceeding the single polaron mo- 
bility at low particle/hole doping and sufficiently strong 
attractive JT interaction, Vjt < — 1 (see Fig. This is 
the most important result of our calculations. The expla- 
nation of this effect follows from the the fact that all par- 
ticles in the cluster are paired, therefore formation of the 
clusters with odd number of particles is highly energeti- 
cally unfavorable. With increased doping and strength of 
the attractive JT interaction also fj,b P starts to show an 
inflection point at the temperature, T M P , which is always 
lower in the bipolaron case. Below (T^ p ) both mo- 
bilities show an Arrhenius temperature dependence (see 
Fig. [2|). The activation energies clearly show a smaller 
energy barrier for bipolaron hopping in comparison to 
single polaron hopping. 

In our simulation ^bp/7o P exceeds ^ S p/7o at small dop- 
ing and low temperature. This however does not im- 
ply that the bipolaron mobility cannot dominate the 
charge transport even at the intermediate dopings and 
for weaker JT attractions. If the binding energy of a 
bipolaron is larger than the phonon frequency, the single 
phonon processes which break the bipolaron are forbid- 
den. The prefactor 70 in the single particle hopping prob- 
ability @ will be therefore strongly suppressed^ because 
more then one phonon is necessary to conserve energy. 
On the other hand, the bipolaron hopping processes do 
not require real bipolaron breaking. As our simulations 
show (see bottom panels in Fig. [5]), the bipolaron hop- 
ping processes can have a lower energy barrier even with 
a moderate JT attraction and can be assisted by a single 
phonon leading to 7q P > 70. This results in a further 
suppression of the single polaron transport with respect 
to the bipolaron transport. 

In the hole doping region at low temperatures the sys- 
tem shows complete orbital (pseudospin) orderJ^ The 
effective attractive JT interaction between holes govern- 



ing the hopping process then becomes isotropic. Since in 
this case the bipolaron mobility also exceed the single- 
polaron one we can conclude that the anisotropy of the 
attractive interaction is unimportant with respect to the 
domination of the bipolaron transport at low T. 



IV. CONCLUSIONS 

Using a new technique for the calculation of the trans- 
port properties of polarons and bipolarons in the charge 
segregated state of doped classical strongly correlated 
systems in the hopping regime, we analyzed the electri- 
cal conductivity of a quasi-2D system with charge segre- 
gated polarons in a textured state driven by the Coulomb 
interaction. The problem of choosing the appropriate 
hopping probability, which has been an issue in previous 
worksi^ has been sidestepped by the fact that the size of 
our bipolaron is small (one lattice constant). We have 
shown that at high temperature the mobility of carriers 
increases with decreasing temperature as T , Single- 
polaron mobility shows a maximum at the characteristic 
temperature where clusters start to form. Below this 
temperature, the single-polaron conductivity is exponen- 
tially suppressed because the single-polaron hopping re- 
quires dissociation of the polaron from a cluster. On 
the other hand, the bipolaron transport shows a weaker 
suppression, which starts at a lower temperature, and 
thus becomes the dominant transport channel at low tem- 
peratures. Apart from resistivity, the Peltier effect, the 
Nernst effect, the Hall effect and particularly the Lorentz 
number are all expected to exhibit a low-temperature 
crossover when bipolaron transport becomes dominant, 
so the model predictions are eminently verifiable by ex- 
periment. 
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